Mean-Field and Perturbation Theory of Vortex-like Composite Fermions 
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We develop a field theory for a partially filled Landau level based on composite fermions with 
a finite vortex core, whose mean-field states are exactly those described by well-tested trial wave 
functions. Despite non-orthogonality of free composite-fermion states and non-Hermiticity of the 
mean-field Hamiltonian, a consistent perturbation theory is formulated and the mean-field Fermi 
sea at half filling is shown to be stable. While the low-energy and long-distance physics is the same 
as in the Chern-Simons fermion theory, new physics is expected to show up for larger wave vectors. 



PACS numbers: 73.40.Hm, 71.10.+X, 71.27.+a, ll.10.Kk. 



Physics of interacting planar electrons in a partially 
filled lowest Landau level (LLL) in a strong magnetic 
field B is surprisingly rich. For example, at the filling 
factor v = 1/0, depending on whether is an odd or 
even integer, the system may be cither in an incompress- 
ible or compressible [Q fluid state, exhibiting the 
quantized |1 or unquantized J(| Hall effect respectively. 
The corresponding ground states are known to be de- 
scribed by numerically well-tested trial wave functions, 
such as the Laughlin and Rezayi-Read Q ones. Phys- 
ically, these wave functions can be understood in terms 
of the composite fermion (CF) scenario ||||: Namely 
at fillings near v = 1/0, with even, a quasiparticle 
is an electron bound to a 0-fold vortex in the electron 
fluid p0| , called a CF, and the strongly correlated fluid 
itself can be described as a collection of such weakly inter- 
acting quasiparticles in a weaker effective magnetic field 
AB = B — B^ with B^ = 0n e 0o. (Here n e is the average 
electron density, 0o the flux quantum.) For example, if 
the residue interactions are ignored, the ground state at 
v = 1/2 corresponds to a Fermi sea of CF's, while that 
at v — 1/3 to a completely filled LLL of CF's in AB. 

In this Letter, we will report on a local field theo- 
retical approach that allows a systematic improvement 
of these trial ground states, which are managed to ap- 
pear as the mean-field ground states. Our approach is 
based on a full realization of the CF scenario, in which 
the CF has a finite vortex core, improving the usual 
Chern-Simons fermion (CSF) theory 11 3,12] in which 
an infinitesimally thin 0-flux is attached to each electron. 
The inclusion of a finite vortex core is implemented by a 
non-unitary transformation, which makes CF states with 
definite momenta non- orthogonal to each other and the 
mean-field and perturbed Hamiltonians non-Hermitian. 
Despite these seemingly troublesome features, it is shown 
possible to formulate a consistent perturbation theory, in 
which the mean-field ground state, a Fermi sea of CF's, 
at v = 1/0 is perturbatively stable. We have also com- 
puted the density-current response functions in the ran- 
dom phase approximation (RPA), and explicitly verified 
that for small wave vectors they indeed agree with those 
in the CSF theory New physics due to the finite vor- 
tex core is expected to show up for larger wave vectors. 



Generalizing a recent composite-boson construction 
]l3] | , we introduce the CF field operators by (with even) 

$(£) = e -'W) ?/,(£), H{x) = ^{x)e J <> (S) , (1) 

where ip{x) is the (spinless) electron field operator, and 

J-Ax) = / d 2 x' p(x') log(z - z') - |z| 2 /4/?, (2) 



with l^ the magnetic length in B^ and z = x + iy. The 
usual CSF transform Q contains only the imaginary part 
of 01og(z— z'), which describes the phases due to a vortex 
with vorticity bound to an electron at z' . We have in- 
cluded the real part too, incorporating a finite vortex core 
and making the transformation non-unitary p4| . Using 

e- J * {S) ij(x') = (z- z'y^{x')e- J i> {3) , 
rj\x')e- J i {S) =(z- zYe- J i {s) yj\x"), (3) 

it is easy to verify {$(£),$(£')} = {n(f),n(f )} = 0, 
and {$(x),IL(x')} = 6~( 2 \x — x'). Obviously <E> and n are 
not Hermitian conjugate: n = & e J * +J ^ . Notice that 
p{x) = ^(x)ip(x) = n(x)$(f) and [J d 2 x'p,n{x)} = 
n(x). Thus the CF density is the same as the electron 
density, and n creates a CF while $ annihilates one. 
In terms of CF, the usual electron Hamiltonian reads 



H = — 



1 



2?7Jfc 



d 2 x n(x)(V +iA- ivi) 2 <f>(x) 



+ - / d 2 xd 2 x' 5p(x)V(x -x')6p(x"), 



(4) 



where m,b is the electron band mass; bp = n$ — n e , 
and v^(x) = «VJ^ = a(x) + in x a(x) — ix/21 2 ^; h is a 

unit vector perpendicular to the plane and a is the usual 
Chern-Simons gauge field, given by 



a(x) = / d 2 x p(x) (n x x)/!^^ 



(5) 



in the gauge V • a = 0, satisfying b — V x a — 27r0p. In 
deriving v^, we have used V(Relogz) = V(Imlogz) x h. 

The physical justification for including the real part of 
01og(z — z') in the CF transformation (0) lies in the fact 
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that the resulting mean-field states give rise to numer- 
ically well-tested wave functions. To see this, we note 
that at the mean field level, one ignores the fluctuations 
of the Chern-Simons field a and, for fillings at or close to 
l/4>, take a to be a classical field determined by eq. (|^) 
with a uniform density n e , as in CSF theory: 



p(x) = n e , a(x) = (B;/2)n xx = Ai(x). 



(0) 



Here one has n x a — which results in J + ,P = 0. 

Substituting eq. @ into eq. (Q), we get the mean- field 
Hamiltonian describing free CF's in an effective field AB: 



Hmf — 



1 



2m;, 



d 2 x II(x)(V + iAA) 2 $(x). 



(7) 



Once we get the CF wave function for a mean-field state, 
x(a?i, xn) = (0|<f>(fi) . . . $(xjv)|MF), the correspond- 
ing electron wave function can be easily read off as 

^mf(x u ...,x n ) = {0\^(x 1 )...^(x N )\MF) 
= (0|e J ^$(:ri) . . . e J ^ N ^{x N )\MF) 

= ~[[( z i - ^)^ ex P[-^2 l Zl ' 2 ] x(xu-,x N ) , (8) 

i<j <fi i 

recovering Jain's rule for trial wave functions . Here we 
used the identity (||) to move all e J ^ Xi ^ to the left, which 
act on the vacuum (0| yielding the Gaussian factor. The 

mean-field CSF theory missed the factor Y\.i<j \ z * ~~ z jl^> 
whose presence in eq. (||) is due to the inclusion of the 
real part of 4>\og(z — z') in our CF transformation (0). 

At exactly v = 1/0, CF's are in zero effective magnetic 
field AB = 0. The mean-field Hamiltonian (Q) requires 
the ground state be a filled Fermi sea of CF's, with Fermi 
wave vector hp — (27r0?i e ) 1 / 2 = 1/^. In terms of the CF 
field in momentum space, the Fermi-sea state ket is 



igo) = n n (£) i°>- 



(9) 



k<k F 



Thus, with xo = {0\<f>(x 1 )...<I>(x N )\G a ) = det(e ifc4 ^), 
eq. (pf) reproduces the (unprojected) Rezayi-Read trial 
wave function 0. According to eq. ([?]), the mean- field 
quantum Hall state at v — p/((f>p + 1) is the state with 
CF's completely filling p Landau levels in AB. For p = 1, 



Xi{xi,-,x N ) = Y[(zi - z^expi-—^— N 2 ] 

i<j , 



(10) 



the mean-field electron wave function (M) just gives the 



Laughlin trial wave function |l|] (since I - + l A 2 B 



Zj y 



exp[ 



4/2 



i- B 2 ) 



(ii) 



To systematically improve the mean- field theory, we do 
perturbation theory to include the effects of fluctuations 



of the Chern-Simons field a. For definiteness, we restrict 
to the filling v = 1/0, for which AA = 0. Instead of eq. 
(0), we consider the following Hermitian Hamiltonian: 



H a + H 1 = 



1 



-1 

2mj 



J d 2 iII(i)[V - i(6a + ifi x 5a)} 2 $(x) 



87T0 2 



d 2 xd 2 x' 5b{x) V(x-x')5b(x"). (12) 



with Ho = Hmf and 5a = a — a. We go from the 
Schrodinger to the interaction picture by a similar but 
non-unitary transformation (Hq is not Hermitian) : 



W)h = e iHot \m)s, Oi(t) = e iHot 6 s e 



-iH t 



(13) 



The CF operators Hi (x, t) and $ / (x, t) satisfy the usual 
canonical equal-time anti-commutation relations. The 
evolution operator, U(t,t') = e iHot e -iH(t-t') e -iH t'^ is 
no longer unitary, but we still have U(t,t) — 1, and 
U(t : t")U{t",t') = U(t,t'), U- l (t,t') = U(t',t). Also the 
Schrodinger equation and the Dyson formula for U(t, t') 
are formally the same as before: 



U(t',t) =T expf-i dTHi(r) 



(14) 



where H\(t) = e lHot Hie~ lHot , and T does time-ordering. 

To proceed, first we need an appropriate basis to eval- 
uate the matrix elements of eq. ([14|) . The base kets and 
bras for CF's with definite momenta ]l5j are given by 

|{£}> = nn$)|o), <{&}| = <o|I[*&)- ( 15 ) 



They are eigenvectors of Hq and have orthonormal over- 
laps. (Note that the bras do not have orthonormal 
overlaps with the kets {fcj}}!) Thus, corresponding to 
the Fermi-sea ket ([)]), the bra describing the CF Fermi 
sea, that has a unit overlap with it, is given by 

(Go | = (0| J] <Go|Go) = 1. (16) 

Both the ket [Go) and the bra (Go \ are eigenvectors of H 
with the same energy e = J2k<k F k 2 /2rrib, while (Go\ is 
not. As a rule, corresponding to usual expectation values 
in the unperturbed ground state, we always consider the 
matrix elements between (Gq\ and \Gq). 
The free CF propogator is defined as 

G (x,t;x',t') = -i(G \T($ I (x,t)Il I (x',t'))\G ), (17) 
whose Fourier transform is the same as a free electron 
9{k - k F ) 6{k F - Jfc) 



u> — tj: — i0 + 



(18) 



where 0(k) is the step function and = k 2 /2mft. In- 
troducing ao to implement the field-density constraint, a 
Lagrangian approach for the Chern-Simons propogator 
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D%(x, t; 0, 0) = -i(G \T{5a^x, t)5a v (Q, 0))|G ) ; (19) 
leads to the one in usual CSF theory 

( 



v{q) 
2ma 



\ 



2ni<f> \ 

q 
o 



(20) 



Here we have adopted the 2x2 matrix formalism Q with 
fx, v = 0, 1; stands for the time component, 1 the space 
component transverse to q. Eq. ( |l2|) implies a complex 
CF-Chern-Simons coupling p5a^ + (J — in x j) ■ 5a, with 
j = (-i/2m 6 )[(VII)$-nV$], resulting in the CF-CF-5a 
vertex (Fig. la) 



9p,(k + q,k) = 1, 



[2k + q) x g + jj(2fc + g) -q 
2m b 



(21) 



where q = q/q is a unit vector. Note that the second 
term in the a = 1 component is absent in the CSF theory. 
Moreover, in view of (5a + in x da) 2 = 0, there is no CF- 
CF-Sa-Sa vertex (Fig. lb) in our Feynman rules. This 
makes the structure of Feynman diagrams more like a 
theory with two-body potential than usual gauge theory. 

Using the Wick theorem and the above Feynman rules, 
one can calculate as usual the matrix elements of the 
evolution operator (|l4|) between free CF states. To make 
connection to physics, we have managed to prove a gener- 
alized Gell-Mann-Low theorem, which relates these ma- 
trix elements to correlation functions: 



(G\T0 1M (x u t 1 )---O n ,H(x n ,t n )\G) 
(G |rn Oij^uUpjoo, -oo)|Gq) 
(G 1 C/(oo, -co) | Go) 



(22) 



where |G) is the true ground state, Oi t n(xi,ti) are lo- 
cal operators in terms of and Hh in the Heisenberg 
picture. The result, central to our paper, implies a consis- 
tent perturbation theory with usual diagrammatic tech- 
niques, despite nonunitarity of our CF transformation 
(|l|). The proof follows the same steps as in usual many- 
body perturbation theory p6[ | , but we have to be careful 
about problems due to non-Hermiticity of Hq and Hi. 

The basis of the theorem (22) lies in the following 
lemma: The state obtained from the free CF Fermi-sea 
| Go) by adiabatically switching on H\(t), 



\G) = C Km ft(0.-°°)!Co) 
^o+(Go|C^(0,-oo)|Go) 



(23) 



is an cigenstate of Hq + Hi which, by the adiabatic hy- 
pothesis, is assumed to be the true ground state of the 
system. Here G is a normalization constant and U^it^t') 
the operator @ with H x (t) -» e'^H^r). 

To make sense, the limit in the right side of eq. ( pjjj ) 
has to exist. We first note that, by purely combinatoric 
considerations as usual, U v = U V Lexp(U v o c ), where U v l 



is the linked part of U v , while U v q c the sum of the con- 
tributions to (Go | U n | Go) from all unlinked connected di- 
agrams. It can be shown that as n — > 0+, U v l(Q, —oo) 
is regular, while [/^(O, — oo) diverges as 1 / 77 due to the 
integration of e vti in eq. (|lj) : 



UnOc(0, -00) = iA/rj + InG, 



(24) 



where A is 77- independent. The potential problem lies 
in the possibility that A may have a non-zero imaginary 
part; then the state U v (0, — oo)|Go) would have either 
zero or divergent norm as 77 — > 0+, making the limit in 
eq. (|23| ) nonsense. Usually U v is unitary, so A is real. 
However, this argument does not apply in our case, since 
our Hi contains a complex Chern-Simons coupling for 
CF's. We have managed to check explicitly up to three 
loops, and to give an argument for any number of loops, 
that diagram by diagram the contribution to A is real. 
Thus eq. (24) gives rise to a harmless divergent phase 
factor to the numerator in eq. (p3|), that is cancelled by 
the denominator, and ensures the perturbative stability 
of the mean-field Fermi-sea ground state |Go). Moreover, 
we have been careful to make sure that the usual proof for 
the state ( p3| ) to be an eigenstate of H + Hi goes through 
in our case: Only the commutation relations are needed 
here; whether Hq and Hi are Hermitian is irrelevant. 
Similarly, by adiabatically switching off Hi(t), we get 



(G| = G lim 



(Go 1 17,(00,0) 



»-o+ (G 1 tT, (00, 0)| Go) ' 



(25) 



and U v (oo,— 00) |Go) = exp(i2A/r))\Go), again ensuring 
the perturbative stability of |G ). With these results, one 
easily arrives at the generalized Gell-Mann-Low theorem. 

In addition to the formal development, we have calcu- 
lated the density-current response functions, K^q^w), 
in the RPA in appropriate limits, to explicitly check that 
our formalism indeed gives reasonable results. The rel- 
evant Feynman diagrams are given in Fig. 2. We note 
that our Feynman rules, compared to usual CSF the- 
ory, have modified the Chern-Simons but not the electro- 
magnetic couplings. Besides, the RPA equation for our 
K^v is different from usual CSF theory. Explicit calcula- 
tions shows that the effects of these two differences can- 
cel for K^v both in the static (lo = 0) and high-frequency 
(lu 3> qkp/mb) limits for small wavevector q <C kp. We 
may summarize the differences between our theory and 
usual CSF theory by rewriting the RPA equation for our 
K^y in the following form: 



K = K°- K°[K° + AK- U- 1 ]- 1 ^ , 



(26) 



where K° is the response function of the noninteracting 
CF system (see Fig. 2), governed by Hq; the matrix AK, 
absent in usual CSF theory, is a diagonal 2x2 matrix 



AK = diag 0, K°n - K° n - (<) 2 /^o°, 



00 



(27) 
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Here K° is the response function of non- interacting CF's 
with two Chern-Simons vertices (Fig. 2); in deriving eq. 
( p6| ) we exploited relations between the explicit expres- 
sions of the one- loop blocks in Fig. 2. Straightforward 
calculation shows that both in the static (ui — 0) and 
high-frequency (w 3> qkp/mi,) limits for small wavevec- 
tor q <C kF, AK just happens to vanish: 



(28) 



Therefore, for long-wavelength fluctuations either in the 
static or high-frequency limit, there is no physical differ- 
ence between our theory and CSF theory for the linear- 
response functions. This provides a consistency check 
of our perturbation theory: On one hand, the size of 
the vortex core in a CF is of the order of the magnetic 
length, so a probe with wavelength much bigger than 
that should not be able to see the core; on the other 
hand, the high-frequency behavior of the response func- 
tion should be constrained by Kohn's theorem JlTj (with 
mass mb), which is indeed satisfied by our RPA results. 
Thus, the experimental predictions for probes with long 
wavelength discussed by Halpcrin, Lee and Read re- 
main unchanged. Furthermore, when higher-order con- 
tributions are included, the "Fermi-liquid corrections" 
mM are needed and can be done as well. The details 



will be given elsewhere 19 1. 

Nevertheless, there is no reason to believe that eq. ( |28| ) 
could be generally true. To uncover the new physics due 
to the finite vortex core in a CF, which is expected to 
show up for shorter wavelengths or larger wave vectors, 
one needs to evaluate the one-loop integrals with param- 
eters in certain intermediate range. Work is in progress. 

To conclude, several comments are in order. First, our 
perturbation formalism can be easily generalized to study 
the incompressible fractional quantum Hall fluids either 
from the CF scenario or from the composite boson sce- 
nario; in the latter the integer (f> in eq. (^|) is taken to 
be odd, resulting a canonical boson pair $ and n [ jl3| . 
Here we have considered an infinite, homogeneous system 
without boundary. It would be interesting to generalize 
the present formulation to a finite compact geometry, say 
a sphere or a torus, which may be helpful for clarifying 
some consequences of the non-orthogonality of CF states, 
such as fractional and mutual exclusion statistics between 
quasiholes and quasielectrons |2^,|2l[ . Also it is worth to 
study the effects due to a finite vortex core for quasipar- 
ticles on the edge of a finite system with boundary. 

Y.S.W. would like to thank F.D.M. Haldane for dis- 
cussions. The work was supported in part by grant NSF 
PHY-9309458. 
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